Hadronization in heavy ion collisions: Recombination and fragmentation of partons 
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We argue that the emission of hadrons with transverse momentum up to about 5 GeV/c in central 
relativistic heavy ion collisions is dominated by recombination, rather than fragmentation of partons. 
This mechanism provides a natural explanation for the observed constant baryon-to-meson ratio of 
about one and the apparent lack of a nuclear suppression of the baryon yield in this momentum 
range. Fragmentation becomes dominant at higher transverse momentum, but the transition point 
is delayed by the energy loss of fast partons in dense matter. 



Data from the Relativistic Heavy Ion Collider (RHIC) 
have shown a strong nuclear suppression of the pion yield 
at transverse momenta larger than 2 GeV/c in central 
An + Au collisions, compared with p + p interactions 0. 
The emission of protons and antiprotons does not appear 
to be similarly suppressed, and the p/ir + ratio reaches 
or even exceeds unity for transverse momenta above 2 
GeV/c These results lack a consistent expla- 

nation in the standard picture of hadron production at 
high transverse momentum, which assumes that hadrons 
are created by the fragmentation of energetic partons. 
Whereas the observed suppression of the pion yield is at- 
tributed to the energy loss of partons during their prop- 
agation through the hot and dense matter created in the 
nuclear collision - a phenomenon commonly referred to 
as jet quenching [|| - the absence of a similar effect in 
the proton spectrum is puzzling 

We propose that hadron production at momenta of 
a few GeV/c in an environment with a high density of 
partons occurs by recombination, rather than fragmenta- 
tion, of partons. Below we show that recombination al- 
ways dominates over fragmentation for an exponentially 
falling parton spectrum, but that fragmentation wins out 
eventually, when the spectrum takes the form of a power 
law. We also show that recombination can explain some 
of the surprising features of the RHIC data, as first sug- 
gested by Voloshin 0. 

In the fragmentation picture Q the single parton spec- 
trum is convoluted with the probability for a parton i to 
hadronize into a hadron h, which carries a fraction z < 1 
of the momentum of the parent parton. It has been ar- 
gued that the fragmentation functions Di^h(z) can be 
altered by the environment |lfj| . The dominant modifi- 
cation mechanism is the energy loss of the propagating 
parton in the surrounding medium, which leads, in first 
approximation, to a rescaling of the variable z. This 
would affect all produced hadrons in the same way, in 
contradiction with the observations at RHIC. 

Another mechanism of hadron production is quark re- 
combination. In the recombination picture, three quarks 



or a quark/antiquark pair in a densely populated phase 
space can form a baryon or meson respectively. The am- 
plitude for this process is determined by the hadron wave 
function. This mechanism has recently been identified as 
the source of unnatural isospin ratios in the production 
of _D-mesons in the fragmentation region in tt~ + A in- 
teractions at Fermilab [llj . Hadron production in heavy 
ion collisions by recombination of quarks has been con- 
sidered before [l2j , primarily at small transverse momen- 
tum. Quark recombination has recently been invoked to 
explain some aspects of the RHIC data, such as the fla- 
vor pattern of elliptic flow 0], and in the context of a 
scaling model |l4| - 

In the formalism of perturbative QCD, recombination 
is a more exclusive process, and falls off faster than frag- 
mentation with increasing transverse momentum. On the 
other hand, for a meson of given momentum P fragmen- 
tation starts out with a parton with much higher momen- 
tum p — P/z, (z < 1), whereas recombination requires 
a quark/antiquark pair where each parton carries only 
about P/2 in average. However, the spectrum of high- 
momentum partons is steeply falling and further reduced 
by the energy loss in dense matter. Our main result is 
that the momentum range, in which recombination pro- 
cesses can successfully compete with fragmentation, may 
extend up to 5 GeV/c in the favorable environment of 
central heavy ion collisions at RHIC. We will show below 
that recombination is effective whenever the phase-space 
distribution of a system of partons has thermal character. 

Let us consider an expanding system of quarks and 
antiquarks. We assume that the recombination of these 
partons into hadrons occurs on a space-like hypersurface 
Sf. The RHIC experiments indicate that the freeze-out is 
very rapid. The measured two-particle correlation func- 
tions are consistent with an extremely short emission 
time in the local rest frame, suggesting a sudden tran- 
sition after which individual hadrons interact only rarely 
|lfil |. In our treatment, we assume that the dense parton 
matter is devoid of dynamical thermal gluons and pre- 
dominantly composed of quarks and antiquarks at the 
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moment of hadronization. We neglect a possible effec- 
tive quark mass, because we are here interested in hadron 
production at high transverse momentum. 

Denoting the density matrix of the parton system on 
Sf as pi , we can express the number of mesons at freeze- 
out as EdN M /d 3 P = (2ir)- 3 (M;P| p { \M\P) . Here 
\M\P) is a meson state with momentum P. We now 



introduce the single-particle Wigner functions for quarks 
and antiquarks, w a {r;p) and w b (r;p), respectively, and 
neglect multi-particle correlations in the density matrix. 
Further introducing the time-like, future oriented unit 
vector u^j r) orthogonal to the freeze-out hypersurface at 
r e Ef [lfj, UJ| > we obtain the following expression for the 
meson emission spectrum: 
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E ^ = lda^^^[dz\^ b (z)\ 2 w a (r;zP+)Mr;(l-^P + )- (1) 
s f a - b o 

Because we are considering the emission of a meson with high momentum, it is most convenient to use a distribution 
amplitude in terms of light-cone coordinates in a frame where P + P~ and the transverse momentum of the meson 
vanishes, z and 1 — z are the momentum fractions carried by the quark and antiquark, respectively, and a, b denote 
the internal quantum numbers of the quarks (spin, flavor, color). For baryons, a similar calculation yields: 
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The overall normalization is fixed by the number of 
quarks on the hypersurface Ef: 



p ■ u(r)w a (r;p) . (3) 



It is important to note that our results are not re- 
ally sensitive to the model used for the recombination 
process. A complete dynamical description in terms of 
QCD is hard to achieve. However, for the observables 
that we discuss below, it is essential to use two com- 
mon features: the probability for the emission of a meson 
(baryon) is proportional to the single parton distribution 
squared (cubed), and the parton momenta sum up to the 
hadron momentum. At this level of sophistication a de- 
scription with the light-cone wave functions replaced by 
three-dimensional, spatial quark wave functions leads to 
identical results. 

To make further progress, we need to specify the par- 
ton Wigner functions. We first consider the case that 
these are described by the exponential tail of a thermal 
distribution w a (r;p) — exp(— (p ■ u — p)/T) with local 
temperature T(r) and chemical potential /i(r), indepen- 
dent of the internal quantum numbers. The constraint 
that the momentum fractions of all partons in the hadron 
wave function must add up to the total momentum P of 
the hadron then insures that the product of all Wigner 
functions entering into the hadron yield is solely depen- 



dent on the hadron momentum: 

w a (r; zP+) w b (r- (1 - z)P+) = exp(-P • u/T), 
w a (r; z x P + ) w b (r; z 2 P + ) w c (r; (1 - Z\ - z 2 )P + ) 

= ex P (-(P.u-/i B )/T), (4) 

where ps = 3/i. One can then perform the integrations 
over the momentum fractions in Q and J2J and obtains 
the result that the baryon-to-meson ratio is independent 
of the momentum P and simply given by the ratio of 
the number of quark degrees of freedom contributing 
to the emission of the specific hadrons: dNs/dNu — 

Z-^a.b.c i Z-^a,b' 

The summations over color, flavor, and helicity give 
rise to degeneracy factors for mesons and baryons, Cm 
and Cb, respectively. In order to derive an upper limit 
to the nucleon-to-pion ratio, we make the simplifying as- 
sumption that all decuplet baryon states A contribute 
to the nucleon yield through decay, but neglect the con- 
tributions to the high-Px pion spectrum from the decay 
of unstable hadrons. If we require the partons to form 
color singlets at recombination, we obtain C w + = 1 and 
C p = 20/(2 x 3!) = 5/3, yielding the result 

dN p /dN„ + = e^' T C p /C^ - | e^/ T . (5) 

Decay of unstable hadron states is likely to reduce this 
value, but will not change the prediction that the p/tt + 
ratio is constant as a function of particle momentum 
and significantly larger than the value expected from 
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quark fragmentation. We note that this result is in good 
agreement with the RHIC data, which show dN p /dN 7T + 
reaching a plateau around unity in the range 2 GeV/c 
< P T < 4 GeV/c 18]. For the A-hyperon yield, we 
need to include the channels A, E° and leading to 
Ca = 4/3. This would predict a plateau at 4/3 in the 
A/if[? yield, which will be diluted by kaons from the de- 
cay of K*. The p/p ratio is also predicted to be indepen- 
dent of momentum and equal to exp(— 2/is/T), again in 
rather good agreement with the RHIC data 

The predictions are radically different, when one con- 
siders a power law spectrum as it is characteristic in per- 
turbative QCD at large transverse momentum: 

w P^(r;p)=A a (l + ^y a . (6) 

One then finds that mesons always dominate 
over baryons at large momentum, dNs /dNM = 
(27 /AP) a C B A 2 JC M Aq, and that eventually parton 
fragmentation wins out over quark recombination. 
For pions, the local ratio of the two contributions 
is dN^z/dN™ c oc P a . On the other hand, for an 
exponential quark spectrum, fragmentation is always 
suppressed with respect to recombination. 

This result constitutes the main insight gained from 
our considerations: Hadron emission from a thermal par- 
ton ensemble is always dominated by parton recombina- 
tion; only when the thermal distribution gives way to 
a perturbative power law at high momentum, does frag- 
mentation become the leading hadronization mechanism. 
The threshold between the two domains depends on the 
size of the emitting system and the hadron species. 

That the meson spectrum from recombination is deter- 
mined, on the average, by w(P/2), whereas the baryon 
spectrum depends on w(P/3), implies that those kine- 
matic properties of the hadron spectrum, which are due 
to collective flow of partons, extend to higher values of 
the transverse momentum for baryons than for mesons 
13|. This effect is clearly visible in the RHIC data 
18L Il9| . which exhibit a linear rise of the elliptic flow 
velocity V2 with Pt, which continues further in Pt for 
protons and hyperons than for pions and kaons. 

Realistic calculations require the specification of the 
freeze-out hypersurface Ef and the parton spectrum. 
Here we assume boost invariance of £f according to 
Tf = z f — const. 01- F° r the thermal part of 

the quark spectrum we use an axially and longitudinally 
expanding thermal source: 

wt{ m y,p T ) = Ath e~^ cosh ^/ T e-y 2 ^ A \ (7) 

characterized by an effective temperature T w 350 MeV, 
which includes the blue shift caused by the radial expan- 
sion, and a rapidity width A w 2. Here ij and y are the 
rapidities in space-time and momentum. The form of 
this spectrum agrees with the results of the parton cas- 
cade VNI/BMS [2(J, which yields a parton distribution 
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FIG. 1: Top: inclusive Pt spectrum of charged hadrons in 
central Au+Au collisions at y / %i = 200 GeV; data taken 
from the PHENIX collaboration. Bottom: Ratio of protons 
to 7r + as a function of Pt- The region below 4 GeV/c is 
dominated by recombination, the region above 6 GeV/c by 
parton fragmentation. 

exhibiting an exponential shape at low transverse mo- 
mentum and a power law shape for high transverse mo- 
mentum. For the power law tail of the parton spectrum 
we choose the results give n by a lowest-order perturba- 
tive QCD calculation [21j, shifted by Apr = —y/\pr 
with A = 1 GeV to account for the energy loss of fast 
partons, and standard fragmentation functions |22|. The 
normalization of the thermal part of the spectrum is ad- 
justed to fit the measured inclusive spectrum of charged 
hadrons from PHENIX [2^ , as shown in the upper frame 
of Fig. ^ The contributions from recombination and 
fragmentation are shown separately to exhibit the loca- 
tion of the rapid crossover between these two mechanisms 
at about 5 GeV/c. We note that in the parton spectrum 
w + w peTt itself, the crossover between the thermal and 
the perturbative part occurs already at about 3 GeV/c, 
consistent with parton cascade predictions |2lJ. The re- 
combination mechanism shifts this point to higher values 
of pt in the hadron spectrum. 

The lower frame of Fig ^ shows our prediction for the 
p/ir + ratio. The rapid drop of its value in the range 
4-5 GeV/c is an unambiguous prediction of our model. 
Experiments at RHIC have not yet been able to probe 
this Pt range, because the identification of protons has 
not been feasible beyond 4 GeV/c. The identification of 
hyperons is possible to higher Pt, and indications of a 
rapid drop in the A/if° ratio have been found 0. 

Figure [21 shows the scaled ratio of particle yields in 
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terized by an exponential momentum spectrum. Such a 
phase may be appropriately called a quark-gluon plasma. 

Note added: We draw attention to the closely re- 
lated recent work of Greco et al. [2(|, who also pro- 
pose that parton recombination can explain the large 
baryon/meson ratio observed at RHIC. 

We thank T. Mehen and M. Asakawa for very help- 
ful discussions. This work was supported by RIKEN, 
Brookhaven National Laboratory, DOE grants DE- 
FG02-96ER40945 and DE-AC02-98CH10886, and by the 
Alexander von Humboldt Foundation. 



FIG. 2: Ratio Raa of fiadron yields in central Au+Au to p+p 
collisions, scaled by the number of binary nucleon-nucleon 
interactions. Data for pions taken from PHENIX |2^l . 



Au + Au and p + p collisions, called Raa- (We use fits 
to data for the p + p yields and our predictions for Au 
+ Au.) The energy loss of fast partons leads to a nu- 
clear suppression in the fragmentation region {Pt > 5 
GeV/c). For low Pt this suppression is counteracted by 
the recombination mechanism, which is absent in p + p 
reactions. Recombination is more important for protons 
than for pions, resulting in much less nuclear suppression 
for protons. 

RHIC data 19j exhibit a strong increase of the 
anisotropic flow parameter w 2 for mesons and baryons 
at small Pt, which finally saturates. This happens ear- 
lier for mesons than for baryons. It has been argued 
that the flow anisotropy originates in the partonic phase 
|25j . In the recombination region mesons at transverse 
momentum Pt reflect the properties of partons with an 
average transverse momentum Pt/2, while baryons re- 
flect those of partons with Pt/3. It follows that vi sat- 
urates later for baryons than for mesons 0. The tran- 
sition to the fragmentation region would provide such a 
mechanism, but it occurs at too high momentum. The 
observed saturation of vi must, therefore, be due to some 
other mechanism or require a more realistic description 
of the space-time evolution of the system. 

In summary, we propose a two component behaviour 
of hadronic observables in heavy ion collisions at RHIC. 
These components include fragmentation of high-pT par- 
tons and recombination from a thermal parton distribu- 
tion. The competition between recombination and frag- 
mentation of partons can explain several of the surprising 
features of the published data. In particular, the proton 
excess at intermediate Pt, the different nuclear suppres- 
sion observed in pion and proton spectra, and the dif- 
ferent saturation thresholds in the elliptic flow, are easily 
explained. We predict that all baryon spectra will exhibit 
a rapid transition around 5 GeV /c to a region dominated 
by parton fragmentation. Finally, our scenario requires 
the assumption of a thermalizcd partonic phase charac- 
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